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ABSTRACT 

This paper presents the characterization of the in-flight beams, the beam window functions and the associated errors for the Planck Low Frequency 
Instrument (LFI). Knowledge of the beam profiles is the key to determining their imprint on the transfer function from the observed to the actual 
sky anisotropy power spectrum. The main beam distortions aff'ect the beam window function, complicating the reconstruction of the anisotropy 
power spectrum at high multipoles, whereas the sidelobes aff'ect the low and intermediate multipoles. The in-flight assessment of the LFI main 
beams relied on the measurements performed during Jupiter observations. By stacking the data from Jupiter transits, the main beam profiles are 
measured down to -20 dB at 30 and 44 GHz, and down to -25 dB at 70 GHz. The main beam solid angles are determined to better than 0.2% at each 
LFI frequency band. To ensure a characterization of the main beam free from the radiometer noise, a dedicated tuning on the Planck pre-launch 
optical model is performed. This approach provides an optical model whose beams fully reproduce the measurements in the main beam region, but 
also allow us to describe the beams at power levels lower than can be reached by the Jupiter measurements themselves. The agreement between 
the simulated beams and the scanning beams is better than 1% at each LFI frequency band. The simulated beams are used for the computation of 
the window functions for the eff'ective beams. The error budget in the window functions was estimated considering both main beam and sidelobe 
contributions, as well as taking into account the radiometer bandshapes. The total uncertainties in the eff'ective beam window functions are: (at 
i ^ 600) 2% and 1.2% at 30 and 44 GHz, respectively; and at ^ ^ 1000, 0.7% at 70 GHz. 

Key words, methods: data analysis - cosmology: cosmic microwave background - instrument: optics 



1. Introduction 

This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration I 2013), de- 



* Corresponding author: M. Sandri, sandri@iasfbo . inaf . it 
^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 



scribes the beams and window functions of the Low Frequency 
Instrument (LFI). 

Detailed knowledge of the instrumental angular response is a 
key requirement for the analysis of high precision measurements 
of the cosmic microwave background (CMB). Current genera- 



countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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tion of experiments employ multi-frequency focal plane arrays 
whose off-axis beams necessarily deviate, to some extent, from 
an ideal, symmetric, Gaussian shape. The radiation patterns of 
every individual detector and their projected angular locations 
need to be reconstructed with great precision to avoid significant 
systematic eff'ects in the data (Hill et al. 2009; Nolta et al. 2009; 
Huffenberger et al. 2010). 

The Planck optical system is designed to ensure high image 
quality over a wide field of view, for detectors spanning over 
1.5 decades in wavelength (Tauber et al. 2010). The LFI opti- 
cal layout is composed of an array of 1 1 corrugated feed horns, 
each coupled to an orthomode transducer which divides the in- 
coming electromagnetic wave into two orthogonal, linearly po- 
larized components. Thus LFI is observing the sky with 11 pairs 
of beams associated with the 22 pseudo-correlation radiometers. 
Each beam of the pair is named LFIXXM or LFIXXS for the two 
polarization states (Main Arm and Side Arm of the orthomode 
transducer, respectively). Here XX is the radiometer chain assem- 
bly number, ranging from 18 to 28. The beams from LFI 18 to 
LFI2 3 are in the V-band (nominally from 63 to 77 GHz); we re- 
fer to them as 70 GHz. The beams from LFI24 to LFI26 are 
in the Q-band (from 39.6 to 48.4 GHz); we refer to them as 
44 GHz. The beams LFI 2 7 and LFI 2 8 are in the Ka-band (from 
27 to 33 GHz); we refer to them as 30 GHz. The optimization 
of the LFI optical system leading to the focal plane configura- 
tion used in flight is described in Sandri et al. (2010), while the 
preliminary characterisation of the LFI beams based on the first 
in-flight data are reported in Mennella et al. (201 1) and Zacchei 
etaL (2011). 

The LFI map-making procedure does not take into account 
the beam profile, which is eff'ectively assumed to be a pencil 
beam. To correct for the beam shape, the angular power spec- 
trum computed from the observed map is divided by the beam 
window function to reveal the intrinsic angular power spectrum 
of the sky. For this reason, beam knowledge directly aff'ects the 
cosmological analyses. Typically, the beam should be mapped 
to less than -30 dB of the peak to achieve 1% accuracy on the 
angular power spectrum (Page et al. 2003). By stacking the data 
from the first four Jupiter transits, the LFI beams have been mea- 
sured down to -20 dB at 30 and 44 GHz, and down to -25 dB at 
70 GHz with an uncertainty of about 0.3% on the angular reso- 
lution and about 0.5% on the main beam ellipticity. In order to 
achieve the beam knowledge at lower power levels and improve 
the accuracy on the angular power spectrum, a substantial eff'ort 
has been made to tune the Planck optical model, presented in 
Tauber et al. (2010), to fit the in-flight measurements of the LFI 
beams. This ensures a good representation of the LFI optics, both 
in terms of main beam and sidelobes. The separation of the in- 
strumental angular response into main beam and sidelobes can 
be somewhat arbitrary. In the framework of this paper, we con- 
sider three regions defined with respect to the beam line-of-sight 
and reported in Fig. 1 : 

1. the main beam, which is defined as extending to 1.9, 1.3, and 
0.9° at 30, 44, and 70 GHz, respectively; 

2. the near sidelobes, which are defined as extending between 
the main beam angular limit and 5°; 

3. the far sidelobes, which are defined as the beam response 
more than 5° from the line-of-sight. 

Almost all the power falls into the main beam region (more 
than 99%). The collected power coming from the region out- 
side the main beam is called straylight and it is one of the ma- 
jor sources of systematic eff'ects in Planck observations, and in 
CMB experiments in general. Straylight impacts the measured 
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Deviation from tine beam axis [degrees] 

Fig.l. Typical shape of a 30 GHz beam (LFI27M). The plot 
shows the distinction between the main beam, near sidelobes and 
far sidelobes. The distinction between "near" and "far" sidelobes 
is of course arbitrary: we mark their boundary at 5°. The peak of 
the spillover of the primary mirror is clearly visible, at an angle 
of roughly 90°. 



signal essentially in two ways: (i) through direct contamination; 
and (ii) in the photometric calibration of the radiometer detected 
signal. The predicted straylight contamination has been detected 
in the LFI maps, and is reported in the companion paper Planck 
Collaboration III (2013). We would like to emphasize that, since 
no direct measurement of LFI sidelobes was performed in-flight, 
an accurate knowledge of the main beams provides a crucial 
means, though indirect, to quantify the straylight contamination, 
as the sidelobes can be estimated by fitting the electromagnetic 
model to the main beam data. 

For the clarity of the present paper and for consistency with 
the Planck companion papers, we need to make clear three im- 
portant definitions: optical beams; scanning beams; and eff'ec- 
tive beams. The optical beam is the optical response of the feed 
horn coupled to the telescope. It is independent both from the 
radiometer response (bandshape and non-linearity) and from the 
satellite motion (spinning and scanning strategy). It represents 
the pure optical transfer function. The main beam properties of 
the optical beams can be evaluated using optical simulations 
performed with methods largely validated by ground measure- 
ments. The scanning beam is the beam that can be directly mea- 
sured in-flight using planet observations. It derives from the op- 
tical beam coupled with the radiometer response, and smeared 
by the satellite motion. So, with respect to the optical beams, 
the scanning beams have slightly larger angular resolution, el- 
lipticity, and solid angle. The effective beam is a beam defined 
in the map-domain, and is obtained by averaging the scanning 
beams pointing at a given pixel of the sky map, by taking into 
account the scanning strategy and the orientation of the beams 
themselves when they point along the direction to that pixel. 
Therefore, whereas for each radiometer there is one correspond- 
ing optical and scanning beam, the same radiometer has as many 
eff'ective beams as the pixels of the observed sky map. The im- 
portance of the eff'ective beams is two-fold: they are used in the 
window function computation; and the solid angles of the ef- 
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Table 1. Approximate dates of the Jupiter observations. The pe- 
riods include the scan by the entire LFI field of view. 



Jupiter transit 


Date 


OD 


Scan 1 (Jl) 
Scan 2 (J2) 
Scan 3 (J3) 
Scan 4 (J4) 


21/10/09-05/11/09 
27/06/10-12/07/10 
03/12/10-18/12/10 
30/07/11-08/08/11 


161 - 176 
410 - 425 
569 - 584 
808-817 



fective beams must be considered in the estimation of the flux 
density of point sources. 

The data analysis pipeline, which starts from Jupiter obser- 
vations and flows down to the window function characterization, 
is discussed in this paper according to the following outline: 
Sect. 2 describes the scanning beams as measured in the first four 
Jupiter transits, and the simulations which provide their best-fit 
model; Sect. 3 describes the effective beams, calculated using 
the simulated beams and taking into account the Planck scan- 
ning strategy; in Sect. 4 we present the LFI window functions. 
An estimate of the propagation of beam uncertainties to the beam 
window functions is reported in Sect. 5. In this section we also 
report the impact of the near and far sidelobes to the window 
function. For the present data release we do not correct the beam 
window function for the sidelobes. Instead their eff'ect is added 
to the total error budget. For the next data release, we plan to in- 
clude a detailed analysis carried out with the in-band integrated 
beams (main beam and sidelobes) that will be included in the 
data reduction pipeline, both in the calibration and in the win- 
dow function estimation. Finally, in Sect. 6 we summarize our 
conclusions. 



2. Scanning Beams 

Jupiter is the best source in the sky for mapping the LFI beams 
with a high signal to noise ratio. The brightness temperature of 
the planet is close to 150 K and gives an antenna temperature 
from 40 to 350 mK depending on frequency, when the dilution 
factor of the beams is accounted for. Because the signal is ex- 
pected to be unpolarized at LFI frequencies, the detector output 
(in antenna temperature Tp,) is proportional to the power func- 
tion of the beam as follows: 

T^io^cpr oc i^(^,0)^p+i^(^,0)^pp+ 

+ {x'''^^ • [m, + \Ei6, 0)^/]} (1) 
TA(e,(f>f oc \E(6,(l>i/^\E(6,(f>i/^ 

+ {;r^^T-[im0)fpP + |^(^,<P]}, (2) 

where E(0,(p)^'^ and E(0,(p)^^ are respectively the co-polar 
and cross-polar electric field components of the beam in the M- 
radiometer and S -radiometer, computed in the main beam frame 
(0 = 9mb and = 0mb); and;^^"^^ is the OMT cross-polarization. 
This parameter was measured during the hardware development 
(D'Arcangelo et al. 2009) and is always less than -25 dB over 
the operational bandwidth, so that the terms between the curly 
brackets are considered negligible. To assess the beam properties 
we used four Jupiter transits named "Jl", "J2", "J3", and "J4". 
Table 1 reports the date and the corresponding observational day 
of each transit. 



2.1. Planet Data Handling 

The LFI in-flight main beam reconstruction is based on a min- 
imisation code described in Burigana et al. (2001) and incorpo- 
rated into the Level 2 Planck LFI DPC pipeline. The code ex- 
ploits the calibrated timelines corresponding to the Jupiter tran- 
sits on the LFI beams, and performs a fit of the beam shape with 
an elliptical Gaussian function. With this Gaussian approxima- 
tion, we have defined the angular resolution in terms of the full 
width half maximum (FWHM), the beam ellipticity (e), and the 
beam orientation (0eii)- Moreover, this fit has been used to de- 
fine the beam center so that the beam pointing directions agree 
with the convention adopted in Planck Collaboration II (2013). 
The fit is performed in the plane of the Planck field of view, 
centred along the nominal line of sight (LOS) defined in Tauber 
et al. (2010). In Fig. 2 the LFI footprint on the sky is reported 
for both polarization arms. The data selection is done using the 
pointing information contained in the satellite Attitude History 
File (Planck Collaboration II 2013), which in turn is used to in- 
fer the nominal LOS direction synchronously with the sampled 
data. The data selected for fits to Jupiter lie on square grids cen- 
tred with respect to the main beam pointing direction, of about 
1.7° in total size at 70 GHz, 2.6° at 44 GHz, and 3.8° at 30 GHz. 

For each radiometer arm, the selected data are characterized 
by an array of data samples specified by the signal amplitude 
(in thermodynamic temperature) during the transit, the positions 
(x, corresponding to the scan circles, and y, corresponding to 
the positions along each scan circle) of Jupiter during the tran- 
sit, and the distances between Planck and the planet itself during 
the transit. An initial guess for the main beam input parame- 
ters and their possible ranges has been evaluated directly on the 
measured timelines, together with an estimate of the noise corre- 
sponding to the sensitivity of the ensemble of signal data. Since 
the noise level has an average value far from zero, an off'set has 
been applied in order to have a noise characterized by a null 
average value. Furthermore, only the data with a signal above 
the 3 cr level from the noise have been considered in the min- 
imisation routine of the fitting code. This implies slightly higher 
error bars, but guarantees a negligible eff'ect due to the back- 
ground. No destriping was performed on the timelines because 
it was found that the l/f noise does not aff'ect the reconstructed 
beam shape above -20 dB. In particular, the l/f noise will not 
aff'ect our later estimate of the window functions because these 
are obtained from simulated beams derived from a tuned opti- 
cal model, as described later in this section. The fit procedure 
gives an analytical description of the LFI beams, through the pa- 
rameters that characterize the elliptical Gaussian profile and the 
corresponding statistical uncertainties; the latter are computed 
using the Minuit processor MINOS^ which calculates the param- 
eter errors by taking into account both parameter correlations 
and non-linearities. 

Table 2 reports the main beam descriptive parameters with 
the estimated uncertainties evaluated from the stacked beams ob- 
tained from the four Jupiter transits. In the bar charts, shown in 
Figs. 3 and 4, the four transits are considered separately and then 
stacked^. It is evident that the four measurements give basically 
the same results. Thus, no optical ageing eff'ects are evident in 
these data, which were taken from October 2009 to August 201 1. 



^ http://seal.web.cern.ch 

' FWHM = ^8xlogio(2)x<,,x<j„; e = crljoi^; (^eii is de- 
fined as the angle between the major axis of the ellipse and the x-axis 
of the LOS frame. 



3 



cd 
o 





LFI 


Beams 


— Main Arms 


0.10 






, 1 1 1 i 1 1 1 1 1 


0.05 


: i 






0.00 




# 21 
% 20 


- 


0.05 




^ ^18 


- 

@ 25 : 


0.10 








-0 


AO 


-0.05 


0.00 0.05 0. 



Planck Collaboration: LFI beams and window functions 

LFI Beams — Side Arms 



cd 
o 



0.10 




1 <<<,,, , 


1 1 1 1 1 1 , 


0.05 


: (9 


" § 22 




0.00 


: i 


# 21 
^ i 19 


- 


0.05 




y ^18 


© 25 : 


0.10 








-0 


.10 


-0.05 0.00 


0.05 0. 



Cross-scan 



Cross — scan 



Fig. 2. Scanning beam profiles for both polarization arms, reconstructed from the first four Jupiter transits. The beams are plotted in 
contours of -3, -10, -20, and -25 dB from the peak at 70 GHz (green), and -3, -10, -20 at 30 (blue) and 44 GHz (pink). 



The improvement in terms of the uncertainties obtained using the 
four scans together is remarkable. 



Table 2. Main beam descriptive parameters of the scanning 
beams, with uncertainties (1 cr). 



14.0 



70 GHz CHANNEL 



Beam 


FWHM 


Ellipticity 


<Aell 


12.0 - 




(arcmin) 




(degrees) 





70 GHz 








18M 


13.41 + 0.03 


1.24 + 0.01 


85.51 +0.68 


18S 


13.47 + 0.03 


1.28 + 0.01 


86.35 + 0.55 


19M 


13.14 + 0.04 


1.25 + 0.01 


78.94 + 0.67 


19S 


13.09 + 0.03 


1.28 + 0.01 


79.12 + 0.58 


20M 


12.84 + 0.03 


1.27 + 0.01 


71.62 + 0.62 


20S 


12.84 + 0.04 


1.29 + 0.01 


72.61 + 0.61 


21M 


12.76 + 0.03 


1.28 + 0.01 


108.00 + 0.52 


21S 


12.87 + 0.03 


1.29 + 0.01 


106.98 + 0.57 


22M 


12.92 + 0.03 


1.27 + 0.01 


102.05 + 0.57 


22S 


12.98 + 0.03 


1.28 + 0.01 


101.74 + 0.57 


23M 


13.33 + 0.03 


1.24 + 0.01 


93.48 + 0.67 


23S 


13.33 + 0.04 


1.28 + 0.01 


93.60 + 0.59 


44 GHz 








24M 


23.23 + 0.07 


1.39 + 0.01 


89.85 + 0.53 


24S 


23.10 + 0.07 


1.34 + 0.01 


89.98 + 0.53 


25M 


30.28 + 0.10 


1.19 + 0.01 


115.41 + 1.02 


25S 


30.92 + 0.10 


1.19 + 0.01 


117.34 + 1.02 


26M 


30.37 + 0.12 


1.20 + 0.01 


62.13 + 1.14 


26S 


30.61 +0.11 


1.19 + 0.01 


61.42 + 1.09 


30 GHz 








27M 


33.06 + 0.10 


1.37 + 0.01 


101.24 + 0.53 


27S 


33.12 + 0.11 


1.38 + 0.01 


101.37 + 0.54 


28M 


33.17 + 0.11 


1.37 + 0.01 


78.53 + 0.57 


28S 


33.28 + 0.10 


1.36 + 0.01 


78.87 + 0.54 




18M 18S 19M 19S 20M 20S 21M 21S 22M 22S 23M 23S 
Radiometer 

30 and 44 GHz CHANNELS 




24M 24S 25M 25S 26M 26S 27M 27S 28M 28S 
Radiometer 

Fig. 3. FWHM at 70 GHz (upper panel) and 30/44 GHz (lower 
panel) for the four Jupiter scans (grey bars) and for the stacked 
beams (white bar). 



2.2. From optical beams to scanning beams 

The optical beams are the optical response of the feed horns cou- 
pled to the Planck telescope. They are independent of both the 
radiometer response (bandshape and non-linearity) and the satel- 
lite motion (spinning and scanning strategy). They represent the 
pure optical transfer function for Planck LFI. 
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1.5 7 




Radiometer 



Fig. 4. Ellipticity at 70 GHz (upper panel) and 30/44 GHz (lower 
panel) for the four Jupiter scans (grey bars) and for the stacked 
beams (white bar). 

2.2.1 . Main beams 

In the main beam region, the optical beams have been evaluated 
from simulations carried out by the application of physical optics 
and the physical theory of diffraction using GRASP^. A dedicated 
optical study has been carried out with the goal of fitting the sim- 
ulated beams to the in-flight measurements. The optical model 
was tuned to minimize the binned residual maps down to -15 dB 
from the power peak, as described in (Planck Collaboration ES 
2013). This approach is preferable to using polynomial fits be- 
cause it is less affected by the noise and the background: the op- 
tical model is much more rigid than polynomial fits, and the full 
focal plane is simultaneously fitted with a single optical model. 
The reason for this study is to produce an ensemble of noise- 
free beams that are representative of the Planck LFI flight optical 
beams, and also both account for beam aberrations at very low 
levels and include the cross-polarization response, which was 
not measured in flight. Of course, before the comparison with 
the data, the optical beams are properly smeared to take into ac- 
count the satellite motion. Beam smearing comes from the fact 
that, while integrating toward a particular direction in the sky, 
the satellite moves and the optical beam is convolved with a top 
hat along the scanning direction. Since during the scanning the 
beam moves rigidly, the convolution is equivalent to an average. 
Whereas this effect is negligible in the calibration step (Planck 
Collaboration V 2013), this is not the case for the main beam 
measurements with planets, for which this effect smears the op- 
tical beam along the scan direction, increasing the beam asym- 
metries in a non-negligible way. 

In Fig. 5 the maps obtained from the difference between 
measurements and simulations for the 70 GHz beams are shown; 
the same comparison has been plotted for the 44 GHz radiome- 
ters in Fig. 6, and for the 30 GHz radiometers in Fig. 7. The 

^ The GRASP software was developed by TICRA (Copenhagen, DK) 
for analysing general reflector antennas (http : //www . ticra . it). 



colour scale is 2.25 times the rms of the beam difference and the 
units of the color bar are in thousandths of the peak height, i.e., 
0.1% of the beam maximum. The colour scale has been sym- 
metrized around in the minimum and maximum values so that 
the zero level is green in all the plots. The size of each patch is 
fixed: 120'xl20' for 30 and 44 GHz and 50'x50' for 70 GHz. 




Fig. 5. Difference between measured (dashed line) scanning 
beams and simulated (solid line) beams (70 GHz channel). The 
colour scale is 2.25 times the rms of the beam difference and the 
units of the colour bar are in thousandths of the peak height, i.e., 
0.1% of the beam maximum. The contours correspond to -3, - 
10, -20, and -25 dB from the peak. The size of each patch is 
50'x50', centred along the beam line of sight. 
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where IImb is the main beam solid angle calculated with the 
pattern normalized to the isotropic level. In the same Table are 
also reported the main beam solid angles of the simulated and 
scanning beams; the values agree to better than 1%. 



Table 3. Beam efficiency and solid angles computed on the op- 
tical beams and simulated beams. In the first column the main 
beam efficiency, t/, derived from optical beams, is reported. In 
the second column we report the percentage of the power enter- 
ing the sidelobes (/si): these values are directly computed as X-t]. 
In the last column the solid angles computed on the scanning 
beams are reported. The 1 cr statistical error on the estimated 
solid angle from the scanning beams is about 0.2%. The com- 
parison between the simulated beams and the scanning beams 
shows that most of the solid angles agree to better than 1%. The 
averaged values of the simulated beams are 194, 850, and 1183 
arcmin^ at 70, 44, and 30 GHz, respectively. The averaged values 
of the measured scanning beams are 193, 849, and 1182 arcmin^ 
at 70, 44, and 30 GHz, respectively. 



Fig. 6. Diff'erence between measured (dashed line) scanning 
beams and simulated (solid line) beams (44 GHz channel). The 
colour scale is 2.25 times the rms of the beam diff'erence and the 
units of the colour bar are in thousandths of the peak height, i.e., 
0.1% of the beam maximum. The contours correspond to -3, - 
10, -20, and -25 dB from the peak. The size of each patch is 
120'xl20', centred along the beam line of sight. 
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Fig. 7. Diff'erence between measured (dashed line) scanning 
beams and simulated (solid line) beams (30 GHz channel). The 
colour scale is 2.25 times the rms of the beam diff'erence and the 
units of the colour bar are in thousandths of the peak height, i.e., 
0.1% of the beam maximum. The contours correspond to -3, - 
10, -20, and -25 dB from the peak. The size of each patch is 
120'xl20', centred along the beam line of sight. 



Table 3 reports the main beam efficiency of each LFI optical 
beam, together with the solid angles. The main beam efficiency 
is defined as: 



2.2.2. Beam validation through deconvolution 

To test the goodness of the beam representation, the maps for 
each individual horn at 30 GHz and 44 GHz have been decon- 
volved using the ArtDeco beam deconvolution algorithm de- 
scribed in Keihanen & Reinecke (2012). The code takes as input 
the time-ordered data stream, along with pointing information 
and the harmonic representation of the simulated beam, to con- 
struct the harmonic Qsim coefficients that represent the sky signal. 
From the harmonic coefficients we further construct a sky map. 
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which is now free from the effects of beam asymmetry, assuming 
that our beam representation is correct. 

Before deconvolution we ran the time-ordered data through 
the Madam map-making code (Keihanen et al. 2010), to remove 
low-frequency noise. We saved the baselines that represent the 
correlated noise component, and subtracted them from the orig- 
inal data stream. The cleaned data thus consist of signal with a 
residual noise component that is dominated by white noise. This 
was then used as input to the deconvolution code. 



LFI26 S1-S2 bmned 



LFI26 S1-S2 bmned 




Fig. 8. Comparison between survey differences of binned and 
deconvolved maps obtained using the simulated beam, LFI 2 6. 
The maps were smoothed to 1° resolution in order to suppress 
noise, and the units are Kelvin. The first row of each image cor- 
responds to the binned map, showing a zoom into the Galactic 
region slightly left from the centre (left), and into an unidenti- 
fied point source at location (-90° ,0°) (right). The improvement 
in the deconvolved images is clear. 

We have run the deconvolution on data from each single sur- 
vey, and looked for residual differences between single- survey 
maps. Results for horn LFI26 are shown in Fig. 8, where the 
difference between first "SI" and second "S2" survey maps is 
reported. The left-hand column shows a zoom into the Galactic 
region at location (-40° ,0°). One image covers a square of width 
13.3°. The right-hand column shows a zoom into a point source 
near location (-90°, 0°). The width of this image is 16.7°. In the 
absence of beam asymmetry and other systematics, the differ- 
ence should be due to noise only. 

The top row shows, for comparison, the difference between 
binned maps. In this case, the maps were binned directly from 
the time-ordered data, without attempting to correct for beam 
effects. A given region on the sky is scanned with different beam 
orientations during the different surveys. This gives rise to the 
residual signal that is evident in the top row images. 

The maps were smoothed to a 1° (FWHM) resolution, in 
order to suppress noise. In the case of binned maps this was 
achieved by smoothing with a symmetric Gaussian beam with 



FWHM of 50'. Combined with the width of the radiometer 
beam, this gives a total smoothing of approximately 1°. 

The bottom row shows the corresponding difference of the 
deconvolved maps. We show the same regions as in the top row 
and with the same scaling. We smoothed the deconvolved har- 
monic coefficients with a 1° (FWHM) Gaussian beam, and con- 
structed a sky map through harmonic expansion. Deconvolution 
almost completely removes the Galactic residual, as well as the 
"butterfly" pattern of the point source. This indicates that the 
simulated beams, based on the tuned optical model, are a good 
representation of the true beams. 

The deconvolution is not part of the nominal pipeline but this 
test provides an important cross-check on the beam representa- 
tion since it tests the beam model against the data in a way that 
is independent from the construction of the model. 

2.2.3. Spectral dependence of beam geometry 

Throughout this work, we have assumed a monochromatic re- 
sponse at each LFI frequency. In fact, the bandpasses are wide, 
and vary in detail from one radiometer to another, even within 
the same band. The effective centre frequency for each band used 
in this paper was calculated assuming a thermal (CMB) spec- 
trum. For different source spectra, the central frequency shifts. 
We must also take account of the fact that the beam pattern has 
some frequency dependence. The geometry of the beams is char- 
acterized by three parameters described in the previous section: 
the FWHM; the ellipticity; and the orientation of the beam j/^eii- 
We have investigated the effect on the LFI beams of assuming a 
power-law spectrum with power index a ranging from -6 to -1-6, 
where a = 2 is representative of the CMB spectrum and a = 
of a flat spectrum. We start by generating GRASP models of the 
main beam By(n) at a set of frequencies defined by splitting the 
bandpass into 15 equally spaced steps centred on the nominal 
central band frequency. The 15 beam maps were then averaged 
by weighting each GRASP map pixel by the bandpass T(y) and 
the source spectrum y^, giving the effective beam pattern 



(4) 



with N = J dv T(y)y^. We then derived the geometric beam pa- 
rameters as a function of a. Since the telescope is achromatic, 
only slight variations of the geometric beam parameters are ob- 
served. In addition, the bandpass averaging process further re- 
duces the variability with respect to the monochromatic case. 
The most interesting result is that the three geometrical parame- 
ters vary nearly linearly with a, with different slopes for each ra- 
diometer. The most sensitive radiometer in FWHM is LFI28-S, 
which had d¥WHM/da of about to 3 x 10""^ degrees. Changing a 
from -h2 to -2 will cause a relative change of at most 0.2%, 0.3% 
and 0.4% respectively in the FWHM at 30, 44, and 70 GHz, weU 
below the 1% level. A similar range of relative variations occurs 
for the beam ellipticity. For the orientation parameter, j/^eih the 
amplitude of di//Q\\/da varies from a minimum value of -2 x 10""^ 
degrees (for the detector LFI24-S) to a maximum of 0.36 de- 
grees (for LFI26-M), so that a change in the spectral index a 
from -2 to -1-2 will produce a rotation of the beam of 1 .4 degrees 
at most, in one direction or the other. Those values (assuming a 
CMB spectrum) contribute to the overall calibration uncertainty 
(Planck Collaboration V 2013), and we emphasise that these un- 
certainties in the beam properties are completely independent 
of the colour corrections needed to adjust intensity scales for 
sources with non-thermal spectra. 
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2.2.4. Sidelobes 

The response of the beam pattern outside the main beam needs 
to be carefully understood, as it may have significant impact on 
the Planck data analysis. Although a full physical optics compu- 
tation could be developed to predict accurately the antenna pat- 
tern of the telescope, this is not the case for the whole- spacecraft 
simulations since the physical optics approach cannot be eas- 
ily applied when multiple diffractions and reflections between 
scattering surfaces are involved. For this reason, we have cal- 
culated the sidelobe patterns through the GRASP multi-reflector 
geometrical theory of diff'raction (MrGTD), which computes the 
scattered field from the reflectors by performing backward ray 
tracing. This represents a suitable method for predicting the full- 
sky radiation pattern of complex mm-wavelength optical sys- 
tems in a reasonable time. The MrGTD sequentially computes 
the diffraction fields from any reffector surfaces that are illumi- 
nated, starting from the feed horn. The sequence of scatterers 
and the type of interaction (reffection or diffraction, occurring 
on each scatter) must be defined in the input to the simulation. 
The simplest (first order) optical contributions producing sig- 
nificant power levels are reflections onto the sub-reflector, onto 
the main reflector, and onto the baflfle, as well as diflfractions by 
the sub-reflector, by the main reflector, and by the baffle. Other 
non-negligible contributions derive from two interactions with 
the reffectors (second order - for example, rays reffected on the 
sub-reflector and then diff'racted by the main reflector), three in- 
teractions (third order - for example, rays reflected on the sub- 
reflector, difl'racted by the main reflector, and then difl'racted by 
the bafliie) and so on. Although MrGTD is, in general, less time 
consuming than a full physical optics calculation, it should be 
applied in a rigorous way in order to obtain reliable results, espe- 
cially at low power levels (down to -50 dBi). In addition, when 
many scattering surfaces are involved, the number of ray tracing 
needed may lead to unacceptable computational time, even with 
MrGTD. Since our analysis requires the production of band- 
integrated patterns to account for the frequency-dependent beam 
responses and the radiometer bandpasses, for now the sidelobes 
simulations have been carried out only up to the first order plus 
two contributions at the second order (reflections and diff'rac- 
tions on the sub-reflector, and then diflfracted by the main reflec- 
tor). 

We found that the beam solid angle in the simulated sidelobe 
region is about a factor two lower than the expected value (see 
/si reported in Table 3). It is evident that the first order approx- 
imation adopted underestimates the actual integrated power in 
the sidelobes. Consequently, the sidelobes model must take into 
account the impact of higher-order contributions combined with 
physical optics analyses. We will do so in the future. 

Careful analysis of the LFI 30 GHz data has revealed the im- 
print of Galactic radiation received through the far sidelobes. 
Such a detection is amplified when taking the diflference between 
maps of even and odd surveys, because, due to the diff'erent satel- 
lite orientation, the coupling between the sidelobe pattern and 
Galactic radiation is reversed. A detailed discussion of the sys- 
tematic eff'ects introduced by sidelobe pickup at 30 GHz is given 
in Planck Collaboration III (2013). The expectations of the side- 
lobe pick-up based on the known level of Galactic emission (as 
measured by Planck itself) and our sidelobe model, are in good 
agreement in morphology with the observed eff'ect, as shown in 
Fig. 9. The residual ring clearly visible in the third panel of Fig. 9 
is in agreement with the previous conclusion that the sidelobe 
model should be improved by adding higher order contributions, 
possibly combined with full physical optics analyses. In the bot- 



tom panel of the Fig. 9 we show the diff'erence between data and 
the simulations, these being amplified by a factor equal to the ra- 
tio of the power entering the sidelobes (computed from the main 
beam efl&ciency) and the integral of the simulated sidelobes. It is 
evident that, once we have re-normalised the sidelobe amplitude, 
the ring artifact almost completely disappears. 

It should be noted that, while the sidelobe eff'ect introduces 
additional complication in the analysis, its detection at 30 GHz 
provides an important validation of the simulated beams, which 
can be trusted even to very low power levels in the higher- 
frequency Planck cosmological channels as well, for which the 
sidelobes signatures are not measurable. This has important con- 
sequences for our future analysis aimed at CMB polarization. 
Our plan is to produce refined in-band integrated beams for each 
radiometer, to be included into the polarization data analysis 
pipeline. 

3. Effective Beams 

The eff'ective beam is the average of all scanning beams that 
cross a given pixel of the sky map, given Plancks scan strat- 
egy. The eff'ective beams capture the complete information about 
the diff'erence between the true and observed images of the sky. 
They are, by definition, the objects whose convolution with the 
true CMB sky produce the observed sky map. Similarly, the ef- 
fective beam window functions capture the diff'erence between 
the true and observed angular power spectra of the sky. We com- 
pute the eff'ective beam at each sky pixel for each LFI frequency 
scanning beam and scan history using the FEBeCoP (Mitra et al. 
201 1) method, as was done in Planck'^ early release (Planck HFI 
Core Team 2011b). 

The pre-computation of the eflfective beams was executed 
at the National Energy Research Scientific Computing Center 
(NERSC) in Berkeley (California). The beam data were deliv- 
ered to the Planck data processing centres (Planck Collaboration 
II 2013; Planck CoUaboration VI 2013) over the network, 
on tape and disk, and ingested into the Data Management 
Component (DMC). FEBeCoP associated application software 
was developed and installed to use the eff'ective beams, e.g., fast 
Monte Carlo full sky convolution codes. 

In estimating the eff'ective beams, a cut-off is applied to the 
input simulated beams. From several tests performed we con- 
verged to a cut-off' radius of 2.5 x FWHM. The beam efficiency 
of the simulated beams within this cut-off' radius is reported in 
Table 4. 

For a detailed account of the algebra involving the eff'ective 
beams for temperature and polarization see Mitra et al. (2011). 
Here we summarize the main results. 

The observed temperature sky J is a convolution of the true 
sky T and the eflfective beam 

f = AQB T, (5) 

where the eff'ective beam can be written for the tempera- 
ture in terms of the pointing matrix A^^ and the scanning beam 
b(rj,p^) as 

Here t represents time samples, A^/ is 1 if the pointing direc- 
tion falls in pixel number /, else it is 0; pt represents the exact 
pointing direction (not approximated by the pixel centre loca- 
tion), and rj is the centre of pixel number j, where the scan- 
ning beam b(rj, p^) is being evaluated (if the pointing direction 



8 



Planck Collaboration: LFI beams and window functions 



Table 4. Main beam efficiencies computed on the simulated 
beams using a cut-olf radius of 2.5 FWHM: main and side arm 
efficiencies are reported in the first two columns, the average 
value of the two arms is reported in the third column, and the 
diff'erence between the two arms is reported in the last colunm. 



Horn 


Main 


Side 


Mean 


(Diff) 


70 GHz 










LFI-18 


0.99345 


0.99262 


0.99304 


-0.00082 


LFI-19 


0.99270 


0.99206 


0.99238 


-0.00065 


LFI-20 


0.99111 


0.99084 


0.99098 


-0.00027 


LFI-21 


0.99115 


0.99105 


0.99110 


-0.00010 


LFI-22 


0.99259 


0.99184 


0.99222 


-0.00075 


LFI-23 


0.99360 


0.99274 


0.99317 


-0.00086 


44 GHz 










LFI-24 


0.99762 


0.99826 


0.99794 


0.00064 


LFI-25 


0.99788 


0.99792 


0.99790 


0.00005 


LFI-26 


0.99787 


0.99793 


0.99790 


0.00006 


30 GHz 










LFI-27 


0.99247 


0.99282 


0.99264 


0.00036 


LFI-28 


0.99230 


0.99284 


0.99257 


0.00054 



falls within the cut-off' radius of 2.5 x FWHM, for LFI chan- 
nels). An analogous formula can be written for the tempera- 
ture + polarization eff'ective beam, including the weight vector 
Wt = [l,7COs(2j^),7sin(2j/r)], as: 



t 



WtW, 



yAtMrj,Pt)WtwJ. (7) 



t 



As an example, in Figs. 11, 12, and 13 we compare images 
of four sources (assumed to be unresolved) from the Planck 
Early Release Compact Source Catalogue (ERCSC) (Planck 
Collaboration VII 2011) and FEBeCoP point spread functions 
(i.e., the transpose of the effective beam matrix) on the same 
patch of the sky for the LFI channel maps. The galactic coordi- 
nates (/, b) of the four sources are shown under the colour bar: 
in our sample, these are, respectively, (305.1°, 57.1°), (86.1°, - 
38.2°), (290.0°, 64.4°) and (184.5°, -5.8°), from left to right in 
the three figures. 

We further performed a 2D Gaussian fit of the effective beam 
at several positions of the sky and studied the distribution of 
the fitted parameters: beam FWHM; ellipticity; solid angle; and 
orientation with respect to the local meridian. In order to per- 
form such statistics, we sampled the sky (fairly sparsely) at 
768 directions which were chosen as HEALpix (Gorski et al. 
2005) A^side=8 pixel centres to uniformly sample the sky. The 
histograms of these quantities are shown in Fig. 14. From the 
histograms, we derive the statistical properties of these quanti- 
ties (mean values and standard deviations), which are provided 
in Table 5). 

In Fig. 10 we show the sky variation of ellipticity, FWHM 
(relative variation with respect to the FWHM of the scanning 
main beam), psi (orientation of the eff'ective beam) and beam 
solid angle (relative variation with respect to the scanning main 
beam solid angle reported in Table 2) of the best-fit Gaussian 
to the eff'ective beam at HEALpix A^side=16 pixel centres for 
70 GHz. The eff'ective beam is less elliptical near the ecliptic 
poles, where more scanning angles synmietrize the beam. 

The main beam solid angle of the eff'ective beam, Q^f^, is 
estimated according to the definition: 4;r 2(^o)/niax(5/j), i.e., 
as an integral over the full extent of the eff'ective beam. From the 



eff'ective beam solid angle we can estimate the eff'ective FWHM, 
under a Gaussian approximation: these are tabulated in Table 5. 
The reported FWHMgff are derived from the solid angles, under a 
Gaussian approximation. The mean(FWHM) are the averages of 
the Gaussian fits to the eff'ective beam maps. The former is best 
used for ffux determination, the latter for source identification. 

Note that the FWHM and ellipticity in Table 5 diff^er sHghtly 
from the values reported in Table 2; this results from the diff'erent 
way in which the Gaussian fit was applied. The scanning beam 
fit was determined by fitting the profile of Jupiter on timelines 
and limiting the fit to the data with a signal above the 3cr level 
from the noise, while the fit of the eff'ective beam was computed 
on the maps of the simulated beams projected in several position 
of the sky; the latter are less aff'ected by the noise. 

In Table 6 we indicate with o!"^^ the beam solid angle esti- 
mated up to a radius equal to the FWHMgff (as defined above), 
while indicates the beam solid angle estimated up to a ra- 
dius equal to twice the eff'ective FWHM (FWHMeff). These were 
estimated according to the procedure followed in the aperture 
photometry code for the Planck Catalog of Compact Sources 
(PCCS) (i.e., if the pixel centre does not lie within the given ra- 
dius it is not included). These additional quantities were eval- 
uated for the production of the PCCS (Planck Collaboration 
XXVIII 2013). 



4. Beam Window Function 

CMB temperature anisotropics are a scalar random field on a 
sphere, and can be decomposed in spherical harmonic coeffi- 
cients: 



where 



r danh)Y*,Jhl Tin) = (^imYtm. 



(8) 



(9) 



The finite angular resolution of an instrument h(h, n!) can be 
described through a convolution in real space: 

To^m^ J dQb(h,n^)T(n^l (10) 

which is equivalent to a low-pass filter in harmonic space, and 
whose eflfective action on the power spectrum can be written as: 



(11) 



where is the beam window function. 

As discussed in the previous section, a basic symmetric 
Gaussian approximation is not a good description of Planck 
eff'ective beams. Moreover, the combination of intrinsic beam 
asymmetry and scanning strategy produces eff'ective beams that 
vary significantly over the sky. Therefore, in order to produce 
accurate estimates of the beam window functions, we are forced 
to use detailed Monte Carlo simulations. This has been imple- 
mented using two approaches: first, full timeline-to-map simula- 
tions, where the CMB signal is convolved with realistic scanning 
beams in harmonic space, and then projected into a TOD through 
the Planck scanning strategy and processed in the same way as 
real data; and second, pixel space convolution of CMB signal 
only maps using the eff'ective beams derived with FEBeCoP. 

In principle, for full- sky maps the eff'ective azimuthally av- 
eraged beam window function can be estimated directly from 
Eq. 11: 

W, = (Cf ^)/Q, (12) 
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Table 5. Mean and standard deviation of FWHM, ellipticity, orientation, and solid angle of the FEBeCoP effective beams computed 
with the simulated beams. FWHMeff is the effective FWHM estimated from the main beam solid angle of the effective beam, ^2eff = 
mean(f2), under a Gaussian approximation. 



Channel 


mean(FWHM) 


o-(FWHM) 


mean(e) 


o-{e) 


mean(</r) 




mean(n) 


o-(n) 


FWHMefF 




(arcmin) 


(arcmin) 






(deg) 


(deg) 


(arcmin- ) 


(arcmin-) 


(arcmin) 


70 


13.252 


0.033 


1.223 


0.026 


0.587 


55.066 


200.742 


1.027 


13.31 


44 


27.005 


0.552 


1.034 


0.033 


0.059 


53.767 


832.946 


31.774 


27.12 


30 


32.239 


0.013 


1.320 


0.031 


-0.304 


55.349 


1189.513 


0.842 


32.34 



Table 6. Band averaged beam solid angles. 



Channel 


(arcmin^) 


spatial var. 
(arcmin^) 


(arcmin ) 


spatial var. 
(arcmin^) 


(arcmin^) 


spatial var. 
(arcmin^) 


70 


200.74 


1.03 


186.26 


2.30 


200.59 


1.03 


44 


832.95 


31.77 


758.68 


29.70 


832.17 


31.81 


30 


1189.51 


0.84 


1116.49 


2.27 


1188.95 


0.85 



where Cj^^ is the power spectrum of simulated CMB-only maps, 
Q is the fiducial model used as input, and the ensemble average 
is taken over the Monte Carlo simulations. However, in a realis- 
tic case we mask out some regions of the sky that are contami- 
nated by foreground, and the above equation no longer applies. 
Instead, using the same notation as in Hivon et al. (2002), we 
can write: 

where the coupling kernel M^' encodes the geometric mode- 
mode coupling effect introduced by masking the sky. However, 
we have verified that for the Galactic mask used for power 
spectrum estimation (Planck Collaboration II 2013; Planck 
Collaboration XV 2013) the differences between full-sky and 
cut- sky window functions are marginal with respect to the er- 
ror envelopes discussed in Sect. 5. Therefore, hereafter we'll use 
full-sky approximation. 

4.1. Timeline-to-map Monte Carlo window functions 

Signal-only timeline-to-map Monte Carlo (MC) simulations are 
produced using Level-S (Reinecke et al. 2006) and HEALpix 
subroutines and the Madam map-maker (Kurki-Suonio et al. 
2009; Keihanen et al. 2010) on the Louhi supercomputer at CSC- 
IT Center for Science in Finland; see Appendix B for details. 

Starting from a fiducial CMB power spectrum we have gen- 
erated a set of sky a^m realizations of this Q that are convolved 
with the beam a^^ obtained from the simulated scanning beams. 
Note that the main beams do not collect the full power of the sig- 
nal, since a small part of the signal spills outside the main beam 
to form sidelobes. In this MC just the main beam up to 4 FWHM 
was simulated, not the sidelobes, so the calculated signal values 
were missing that part of the power that goes to the sidelobes. 
This was taken into account at the map-making stage. Note that 
the main beam definition used here (4 FWHM) differs from that 
adopted in the effective beams computation (2.5 FWHM). The 
consequences are discussed in Sect. 5. 

The CMB timelines for each realization were produced ac- 
cording to the detector pointing for each radiometer, and maps 
were made from these CMB timelines with Madam. The same 
Madam parameter settings were used as for the flight maps 
(Zacchei et al. 2011; Planck CoUaboration II 2013). The caU- 



bration step was not simulated, as the simulated signal was con- 
structed as already calibrated, except for the effect of power lost 
to the sidelobes. The impact of sidelobes on the calibration of 
flight data is discussed in Planck Collaboration V (2013), while 
its effect on the beam window function will be discussed in 
Sect. 5 of this paper. For the MC, we assumed that the calibration 
compensates for the missing power in the main beams, according 
to the discussion in Sect. 2.2 of Planck Collaboration V (2013). 

We produced in this way 30 GHz, 44 GHz, and 70 GHz 
frequency maps, and the "horn-pair" maps for 70 GHz 18/23, 
19/22, and 20/21 from the 15.5 month nominal survey. The com- 
putational cost of producing one realization of this set was about 
2000 CPUh. Given this relatively large computational cost, we 
have generated only 102 CMB realizations. Although this leaves 
some residual scatter in the estimated beam window functions 
especially at low multipoles, these maps have been generated 
mostly as a consistency check with respect to the FEBeCoP ap- 
proach as described below, and therefore the number of simula- 
tions is adequate for this purpose. 

Full-sky, timeline-to-map Monte Carlo based beam window 
functions are shown in Fig. 15 for 30, 44, and 70 GHz frequency 
maps. For 70 GHz we also show the beam window functions ob- 
tained considering only subsets of detectors, namely LFI 18-2 3, 
LFI19-22, and LFI20-21. 



4.2. FEBeCoP window functions 

FEBeCoP beam window functions are shown in Fig. 16 for 30, 
44, and 70 GHz frequency maps. For 70 GHz we also show the 
beam window functions obtained considering only subsets of 
detectors, namely LFI18-23, LFI19-22, and LFI20-21. These 
are computed using the effective beams obtained from the simu- 
lated scanning beams with a cutoff radius of 2.5 FWHM. The 
resulting window functions in full sky approximation are ob- 
tained by averaging Eq. 12 over 1000 signal only simulations, 
where every simulated CMB maps is convolved with the effec- 
tive beams described in Sect. 3. 

In Fig. 17 we show a comparison between MC-based and 
FEBeCoP beam window functions. Although there are some 
high-^ discrepancies at 70 GHz, these are located at ^ > 1300 
where the beam W^s drop below 0.01. 
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In addition, however, as we will explain in the next section, 
we should account for the effect of the different choices for the 
cutoff radius between the two methods. Our FEBeCoP calcula- 
tions used a 2.5 FWHM cutoff radius for the main beam, while 
the timeline-to-map, Monte Carlo window functions are derived 
using a 4 FWHM cutoff. In order to quantify the agreement, we 
also show in Fig. 17 the ±lcr error envelopes that will be dis- 
cussed in the next section. So we can fairly conclude that in the 
"region of interest" the two methods agree to 1% level. Since 
the FEBeCoP algorithm is faster than the timeline-to-map Monte 
Carlo, it allows for a significantly larger number of simulations 
(1000 vs 102), resulting in a more accurate estimation of the 
window functions. For the same reason, FEBeCoP also allows to 
perform a robust error assessment as presented in the next sec- 
tion. Hence, the FEBeCoP window functions will be used for the 
power spectrum analysis (Planck Collaboration II 2013; Planck 
Collaboration XV 2013), and will be distributed within the data 
release. 

With FEBeCoP we also estimate the level of contamination of 
the transfer functions due to a non-uniform sky sampling within 
the pixels, comparing the ideal HEALpix pixel window func- 
tion, which is derived under the assumption of uniform cover- 
age, with the true one computed with FEBeCoP. This effect acts 
as a noise term in the maps and it becomes important only at 
very high An analytic treatment of the contamination on the 
maps is described in Appendix F of the HFI Beam paper (Planck 
Collaboration VII 2013). We quote here the level of distortion 
of the window function: 0.1% at ^ = 600 for 30 GHz, 0.4% at 
£ = 800 for 30 GHz, and 0.5% at ^ = 1400 for 70 GHz, in all the 
cases within the error bars. 

5. Error Budget 

We discuss here the main sources of uncertainties in the window 
functions that have proven to be relevant for the LFI. 

5.1. Main beam knowledge 

The propagation of the uncertainties in the beam knowledge to 
the window function has been carried out with a dedicated MC 
pipeline on the Planck optics. The tuned optical model (Planck 
Collaboration ES 2013) was used as the basis to run MC sim- 
ulations with about 500 realizations of the Planck optics. More 
specifically, the wavefront at the aperture of the telescope has 
been artificially distorted by adding to the primary reflector 
randomly-varying amplitude distortions described as modes of 
Zernike polynomials, up to the fifth order. 

The idea behind this assumption is that the true flight beam 
comes from a true flight field distribution at the telescope aper- 
ture that gives the true wavefront. Any small difference between 
our telescope model and the real one can be mapped by aber- 
rations on the aperture wavefront. For each wavefront we sim- 
ulated using GRASP the corresponding beam and we selected 
only the beams with parameters (angular resolution, elliptic- 
ity, and beam orientation) in line with those measured in flight 
within 3cr level. We used the errors in the determination of 
each parameter reported in Table 2. We repeated this for all the 
twenty-two LFI beams producing a set of 3036 beams (corre- 
sponding to 138 slightly different optical models). In this way 
we built a set of beams calculated from plausible optical models 
of the telescope and whose parameters are in agreement with the 
parameters we measured in flight. Then this set of beams was 
used as input to FEBeCoP to compute the corresponding window 
functions. 



The three parameters used in the comparison between sim- 
ulations and measurements (angular resolution, ellipticity, and 
beam orientation) are strongly correlated and this original 
method to obtain the errors on the window function using a 
MC pipeline on the optics takes this correlation properly into 
account, avoiding unphysical solutions in which no correlation 
is assumed. In Fig. 18 we show the beam window functions at 
70GGz for all the 138 simulated optical models. Window func- 
tions are normalized to the fiducial for the 70 GHz channel. 

5.2. Cutoff radius in tlie main beam computation and impact 
of sidelobes 

The impact of sidelobes on the calibration has been discussed in 
Planck Collaboration V (2013). The main result for the discus- 
sion presented here is that the gain values are unbiased, and this 
imposes a constraint on the dipole term of the window function, 
i.e., Wi = I, which fixes the normalization. In principle, in order 
to fully account for beam effects in the window function, one 
should perform a computation of the window function including 
the full beam pattern, with either FEBeCoP or the timeline-to- 
map Monte Carlo. However, this would have a huge computa- 
tional cost making it unfeasible in a Monte Carlo approach. As 
a result, LFI beam window functions are derived from Monte 
Carlo simulations including the main beam only, and therefore it 
is important to assess the effect of neglecting sidelobes. 

We have done a preliminary evaluation by using a analyti- 
cal approach for calculating the window function. In fact, for a 
given azimuthally synmietric beam profile bsiO), the correspond- 
ing ^- space function can be computed using the Legendre 
transform: 

= = 2nj dcos(6)h(6)Pm, (14) 

and = B^. We have then considered the symmetrized beam 
profile for LFI18M detector, and we have computed the corre- 
sponding B{ cutting the integration at different angles, namely 
2.5 FWHM, 4 FWHM, 2°, and 3°, and imposing the normaUza- 
tion constraint at ^ = 1. We report in Fig. 19 the relative differ- 
ence between B^s obtained for the four cutoff values and the one 
resulting from full integration. As expected, results show a small 
^-dependent correction aflrecting mostly the large angular scales. 
In particular, extending the calculation up to the near sidleobes 
makes this effect negligible. 

As already stated in Sect. 4, the combination of intrinsic 
asymmetries in the Planck beam and scanning strategy forced 
us to discard a simple beam symmetrization. The same argu- 
ment applies here, especially considering that near and far side- 
lobes are even more asymmetric than the main beam. Therefore, 
we have extended FEBeCoP calculation for the same LFI18M de- 
tector to cutoff radii of 4 FWHM and 3° (corresponding to ~ 9 
FWHM). In this case we have used the latter window function 
as a reference to compare with when computing the relative dif- 
ference. Results are reported in Fig. 20, confirming that similar 
conclusions can be drawn for a realistic case as well. 

A further improvement in the assessment of the sidelobes ef- 
fect on the window function has been carried out considering the 
variation across the band of the sidelobes themselves. Whereas 
the impact of the main beam variation across the band is small, 
this is not true for the near and far sidelobes. The 47c beams 
of the radiometer LFI IBM (main beam, near and far sidelobes) 
were computed at about twenty frequencies across the radiome- 
ter bandpass and they were averaged taking into account the ra- 
diometer bandshape. The resulting averaged beam has been used 
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to evaluated the impact on the beam window function using the 
analytical approach described above. The shape of the bias is 
very close to that reported in Fig. 19 but the amplitude is slightly 
different with respect to the monochromatic beam. In the error 
budget we have considered the worst case in order to be conser- 
vative as far as possible. 

5.3. Total error budget on window functions 

Using the set of simulated beam window functions we have also 
built the covariance matrix C in ^-space computing: 

= {(W, - {W,})(Wr - (Wr))} , (15) 

where the averaging is performed on the 138 simulations. Then 
we have decomposed in eigenvalues (A^^) and eigenvectors (V^) 
the covariance matrix. In Fig. 21 we show the first eigenmodes 
for the 70 GHz channel. All the error content is substantially en- 
compassed by the first two eigenvalues, that account for cutoff* 
radius and main beam uncertainties respectively. Figures 22 and 
23 show the eigenmodes for the 44 and 30 GHz, respectively. 
The eigenmodes can be used as input of the Markov Chain Beam 
Randomization (MCBR) marginalization code to account for 
beam errors in cosmological parameter estimation (Rocha et al. 
2010). We apply the MCBR procedure to a simulated 70 GHz 
dataset and we find that the parameters mostly aff'ected are n^, 
Qb^^, ^ch^, and As; the increase of the errors can be quantified 
respectively as 12% of cr for the first and less then 8% of cr for 
the others. 

6. Conclusions 

The optics and electronics of the Planck detectors, combined 
with the satellite motion, determine the instrumental angular re- 
sponse to a sky signal. An accurate characterization and a thor- 
ough understanding of the LFI beam patterns is the key to deter- 
mining their imprint on the transfer function from the observed 
to the true sky anisotropy spectrum. In this paper we discussed 
the algorithms used to reveal the most significant LFI beam fea- 
tures that impact the exploration of the underlying cosmology. 
The in-ffight assessment of the LFI main beams relied mainly 
on the measurements performed during Jupiter crossings. The 
calibrated data from four Jupiter scans were used to determine 
the so called scanning beams: the signal-to-noise ratio for this 
data is such as to make it possible to follow the LFI beams 
profile down to -20 dB from the peak, corresponding to dis- 
tances from the beam line of sight of about 1.25 FWHM, i.e., 
the inner parts of the main beams. Fitting the main beam shapes 
with an elliptical Gaussian, we could express the uncertainties 
of the measured scanning beam in terms of statistical errors for 
the Gaussian parameters: ellipticity; orientation; and FWHM. 
While this method allows the accurate in-flight measurement of 
the LFI main beams, the (lower) angular response at larger dis- 
tances from the beam centroid (near and far sidelobes) cannot 
be directly measured from a single point source signal, mostly 
because of the noise and background donunance, so it must be 
modelled diff'erently. Therefore, a further step has been taken to 
build an optimal model for the full LFI beams profile. We devel- 
oped a tuned optical model such that the simulated beams would 
provide the best fit to the available measurements of the LFI 
main beams from Jupiter: we found that this model represents all 
the LFI beams with an accuracy of 1%, which has been consid- 
ered in the propagation of the uncertainties at the window func- 
tion level. The corresponding simulated sidelobes are, in turn. 



consistent with the eff'ect induced by the Galactic spillover as 
observed in survey maps diff'erences. This model, together with 
the pointing information derived from the focal plane geometry 
reconstruction (Planck Collaboration II 2013), gives the most ad- 
vanced and precise noise-free representation of the LFI beams. 
These were also independently cross-checked through a beam 
deconvolution test. The simulated beams were the input to cal- 
culate the eff'ective beams, which take into account the specific 
scanning strategy to include any smearing and orientation eff'ects 
on the beams themselves. The approach was validated by com- 
paring the eff'ective beam Point Spread Functions with images 
from the Planck Catalogue of Compact Sources. 

To evaluate the beam window function, we adopted two in- 
dependent approaches, both based on Monte Carlo simulations. 
In one case, we convolved a fiducial CMB signal with realis- 
tic scanning beams in the harmonic space to generate the corre- 
sponding timelines and maps; in the other case, we convolved 
the fiducial CMB map with eff'ective beams in the pixel space. 
The two methods agree to 1% level. 

To evaluate the error on the resulting window functions, we 
took into account the fact that they were calculated assuming 
full-power main beams. Thus, part of the error budget comes 
from the propagation of the main beam uncertainties through- 
out the analysis, while another contribution comes from neglect- 
ing near and far sidelobes in the Monte Carlo simulation chain. 
We found that the two error sources have diff'erent relevance de- 
pending on the angular scale. Ignoring the near and far side- 
lobes is the dominant error at low multipoles, while the main 
beam uncertainties dominate the total error budget at ^ > 600. 
Representative values of the total error, for scales of cosmo- 
logical interest, range from 0.3% (£ ^ 200) to about 0.8% 

^ 1200). The total uncertainties in the eff'ective beam win- 
dow functions are: 2% and 1.2% at 30 and 44 GHz, respectively 
(at £ « 600); and 0.7% at 70 GHz at ^ « 1000. 
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Appendix A: LFI beams notation 

In table A. 1 we report a selection of the most important symbols 
used in this paper. 

Appendix B: Timeline-to-map Monte Carlo 
simulations 

Signal-only timeline-to-map Monte Carlo (MC) simulations are 
produced using Level-S (Reinecke et al. 2006) and Healpix 
(Gorski et al. 2005) subroutines, and the Madam map-maker 
(Kurki-Suonio et al. 2009; Keihanen et al. 2010) on the Louhi 
supercomputer at CSC-IT Center for Science in Finland. 

Starting from a fiducial CMB power spectrum with ^max = 
3000, we used syn_alm_cxx to generate a set of sky realiza- 
tions of this Q. Starting from the simulated scanning beams we 
calculated their beam a^m using beam2alm with beam_lmax = 
5400 and beamjnmax = 14. The sky a^m were convolved with the 
beama^^ using conviqt_v3 with conv_lmax = 3000, lmax_out 
= 3000, beammmax = 14. The output is a "ringset" table for each 
realization, i.e., a grid of observed sky signal values for 6001 
values of 0, 3002 values of 6 and for 29 beam orientations (</r). 
Note that the main beams do not collect the full power of the sig- 
nal, since a small part of the signal spills outside the main beam 
to form sidelobes. In this MC just the main beam was simulated, 
not the sidelobes, so the calculated signal values were missing 
that part of the power that goes to the sidelobes. This was taken 
into account at the map-making stage. 

The CMB timelines for each realization were produced with 
multimod. The detector pointing (0, 6, i//) for each radiometer 
was reconstructed internally using satellite pointing information 
and the focal plane geometry. The observed CMB signal for each 
sample was interpolated in and from the "ringset" table for 
its pointing using interpol.order = 9 (the effect of beam ori- 
entation is solved exactly by multimod for the beammmax rep- 
resentation of the beam) and output to CMB timeline files. The 
detector pointing reconstructed by multimod was also output 
to disk for the use of the Madam map-maker. The reconstructed 
pointing was compared to the pointing used at the LFI DPC for 
the flight data maps to check that their agreement was satisfac- 
tory. There is an option in multimod, "sampler", to simulate the 
scanning motion of the radiometer during measurement of one 
observation sample. This, however, increases the computational 
cost so much that we turned this sampler off", and simulated this 
scanning motion by using the scanning beams (see above) in- 
stead of the optical beams. 

Maps were made from these CMB timelines with Madam us- 
ing the reconstructed pointing. The same Madam parameter set- 
tings were used as for the flight maps, see Zacchei et al. (2011) 
and Planck Collaboration II (2013). As already discussed in 
Sect. 4, the calibration step was not simulated. 
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Table A.l. Selected LFI beams analysis notation. 



Symbol Description 



M,S 
LOS frame 

U, V 

e 

FWHM 

^ell 
^ max ' max 

T(y) 
Ei (<) 

-^omt 
Pi 

Ati 
Bij ^ B 
b(r,p) 

y 

w 



'opt 



^ff 

FWHMefF 

We 



Main and Side radiometer arm 

telescope's Line of Sight reference frame 

polar coordinates in the LOS frame 

polar coordinates in the main beam frame 

cartesian dimensionless coordinates in the LOS frame 

beam ellipticity 

full width half maximum 

beam orientation defined with respect to the x-axis of the LOS frame 

polarization angle (angle between the detector's polarization axis and the local meridian) 

standard deviation of the elliptical Gaussian 

Main beam efficiency 

bandpass 

detector output in antenna temperature for the M (S) radiometer in the main beam frame 

co-polar electric field component of the beam in the M (S) radiometer 

cross-polar electric field component of the beam in the M (S) radiometer 

cross-polarization of the orthomode transducer 

pointing direction for time sample t 

/-pixel center direction 

Pointing matrix for pixel / and time sample t 

Eff'ective beam for pointing pixel / and beam pixel j 

Scanning beam at a direction f = [6, 0] with the pointing angles p 

optical beam profile 

polarization efficiency y = (1 - 6)/(l + e), being e the cross-polar leakage 

polarization weight factor 

solid angle of the scanning beam 

solid angle of the optical beam 

solid angle of the simulated beam 

solid angle of the eff'ective beam 

eff'ective beam full width half maximum 

beam window function 

eigenvalues of the covariance matrix 

eigenvectors of the covariance matrix 
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Fig. 9. Survey 1 - Survey 2 difference maps. Top: 30 GHz map, 
in yuK. Second: simulated Galactic straylight. Third: data minus 
simulations. The ring still visible in the third panel suggests that 
the model should be improved by adding higher order contribu- 
tions, possibly combined with physical optics analyses. Bottom: 
difference between data and simulations, amplified by a factor 
equal to the ratio of the power missing the main beam (/si) and 
the power entering the simulated sidelobes using the first order 
approximation (in this case, about 1.93). The grey-band is a lack 
of data. 
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Fig. 10. Main parameters of the LFI effective beams: ellipticity 
(first row); FWHM (relative variation with respect to the FWHM 
of the scanning main beam, second row); il/ (third row); and 
beam solid angle (relative variation with respect to the scanning 
main beam solid angle reported in Table 2, fourth row), for the 
70 GHz channel. 
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Fig. 11. Four ERCSC sources as seen by LFI 70 GHz channel (upper panel); linear scale FEBeCoP PSFs computed using input 
simulated beams (central panel); both in arbitrary units. Bottom panel: PSF iso-contours shown in solid line, elliptical Gaussian fit 
iso-contours shown in broken line. PSFs are shown in log scale. 
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Fig. 12. Four ERCSC sources as seen by LFI 44 GHz channel (upper panel); linear scale FEBeCoP PSFs computed using input 
simulated beams (central panel); both in arbitrary units. Bottom panel: PSF iso-contours shown in solid line, elliptical Gaussian fit 
iso-contours shown in broken line. PSFs are shown in log scale. 
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Fig. 13. Four ERCSC sources as seen by LFI 30 GHz channel (upper panel); linear scale FEBeCoP PSFs computed using input 
simulated beams (central panel); both in arbitrary units. Bottom panel: PSF iso-contours shown in solid line, elliptical Gaussian fit 
iso-contours shown in broken line. PSFs are shown in log scale. 
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Fig. 14. Histograms of the three fit parameters (beam FWHIM, ellipticity, and orientation with respect to the local meridian) for the 
eff'ective beams computed using FEBeCoP with the simulated beams. We sampled the sky (fairly sparsely) at 768 directions, chosen 
as HEALpix Nside=8 pixel centres to uniformly sample the sky. 
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Fig. 15. Timeline-to-map, Monte Carlo-based beam window 
functions for Planck 30, 44, and 70 GHz frequency maps. For 
70 GHz, we also show the beam window functions for a subset 
of paired horns, namely LFI18-23, LFI19-22, and LFI20-21. 
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Fig. 16. FEBeCoP beam window functions for Planck 30, 44, and 
70 GHz frequency maps. For 70 GHz, we also show the beam 
window functions for a subset of detectors, namely LFI 18-2 3, 
LFI19-22, and LFI20-21. 
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Fig. 18. Beam window functions generated and normalized to Fig. 21. First four eigenmodes of the covariance matrix of the 
the fiducial. 70 GHz channel. 
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Fig. 19. Relative diff'erence between B^s computed for various 
cutoff' values and the one with full integration. All the functions 
have been computed using eq. 14 for a symmetrized version of 
the LFI IBM beam profile. 
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Fig. 20. Same as Fig. 19, but the window functions have been 
computed with FEBeCoP. 



Fig. 22. First four eigenmodes of the covariance matrix of the 
44 GHz channel. 
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Fig. 23. First four eigenmodes of the covariance matrix of the 
30 GHz channel. 
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